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We investigate low energy nucleon dynamics in the effective field theory (EFT) of nuclear forces. In 
leading order of the two-nucleon EFT we show that nucleon dynamics is governed by the generalized 
dynamical equation with a nonlocal-in-time interaction operator. This equation is shown to open 
new possibilities for applying the EFT approach to the description of low energy nucleon dynamics. 
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I. INTRODUCTION 



Understanding how nuclear forces emerge from the fun- 
damental theory of quantum chromodynamics (QCD) is 
one of the most important problems of quantum physics. 
To study hadron dynamics at scales where QCD is 
strongly coupled, it is useful to employ effective field the- 
ories (EFT's) [1,2], invaluable tools for computing physi- 
cal quantities in the theories with disparate energy scales. 
Following the early work of Weinberg and others [3-7], 
the EFT approach has become very popular in nuclear 
physics [8,9]. A fundamental difficulty in an EFT descrip- 
tion of nuclear forces is that effective Lagrangians which 
are used within its framework yield graphs which are di- 
vergent, and give rise to singular quantum-mechanical 
potentials. These potentials do not satisfy the require- 
ments of ordinary quantum mechanics and need to be 
regulated, and renormalization must be performed. In 
this way one can successfully perform calculations of 
many quantities in nuclear physics. However, in this case 
one cannot parametrize the interactions of nucleons, by 
using some Lagrangian or Hamiltonian. In fact, in quan- 
tum field theory, knowing the Lagrangian is not sufficient 
to compute results for physical quantities. In addition, 
one needs to specify a way to make all infinite quantities 
finite. Another consequence of this situation is that there 
are not any equations for renormalized amplitudes in sub- 
tractive EFT's: In order to resolve this problem, in some 
EFT's [6,10,11] finite cutoff regularization is used. How- 
ever, renormalization is required to render such theories 
consistent, and certain cutoff-dependent terms have to 
be absorbed into the constants before determining them 
from empirical data. 

If the EFT approach, as it is widely believed, is able 
to provide a fundamental description of the interactions 
of nucleons at low energies, one can hope that it will give 
rise to parametrization of these interactions by interac- 
tion operators as fundamental as the Coulomb potential, 
that parametrizes the interaction of charged particles in 
low energy QED. In the quantum mechanics of particles 
interacting via the Coulomb potential, which is an exam- 
ple of the effective theory, one deals with a well-defined 
interaction Hamiltonian and the Schrodinger equation 
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governing the dynamics of the theory. This theory is 
self-consistent, and provides an excellent description of 
atomic phenomena at low energies. In light of this fact 
one can expect that the EFT of nuclear forces will allow 
one to construct a well-defined operator parametrizing 
the interaction of nucleons and governing their dynam- 
ics. At the same time, since the EFT leads to singular 
nucleon potentials, this operator must not be an inter- 
action Hamiltonian. In the present paper we show that 
recent developments in quantum theory provide the pos- 
sibility of a consistent description of nucleon dynamics 
predicted by the EFT approach. 

The above problem of subtractive EFT's is the same 
that arises in any quantum field theory with UV diver- 
gences: Regularization and renormalization allow one to 
render the physical predictions finite, however, it is im- 
possible to construct a renormalized Hamiltonian acting 
on the Fock space, i.e., after renormalization the dynam- 
ics of the theory is not governed by the Schrodinger equa- 
tion. This equation is local in time, and the interaction 
Hamiltonian describes an instantaneous interaction. On 
the other hand, locality is the main cause of UV diver- 
gences in quantum field theory (QFT), and hence regular- 
ization and renormalization may be considered as some 
ways of nonlocalization of the theory. In Ref. [12] it 
has been shown that the use of the Feynman approach 
[13] to quantum theory in combination with the canon- 
ical approach allows one to extend quantum dynamics 
to describe the evolution of a system whose dynamics 
is generated by a nonlocal-in-time interaction, and an 
equation of motion has been derived as the most general 
dynamical equation consistent with the current concepts 
of quantum theory. Being equivalent to the Schrodinger 
equation in the particular case where interaction is in- 
stantaneous, this equation permits the generalization to 
the case where the interaction operator is nonlocal in 
time. It has been shown [12] that a generalized quantum 
dynamics (GQD) developed in this way provides a new 
insight into the problem of UV divergences. 

The aim of the present paper is to show that the for- 
malism of the GQD developed in Ref. [12] opens new 
possibilities in the effective theory of nuclear forces. We 
show that in leading order of the two-nucleon EFT nu- 
cleon dynamics is governed by the generalized dynami- 
cal equation with a nonlocal-in-time interaction opera- 
tor. By using the example of the T matrix that describes 
the contact term in leading order, we investigate the dy- 
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namical situation in an EFT after renormalization. We 
show that this T matrix has the properties that are com- 
pletely unsatisfactory from the point of view of ordinary 
quantum mechanics, but satisfies the generalized dynam- 
ical equation. Moreover, in this case we deal with the 
dynamics that is described by the model developed in 
Refs. [12,14] as a test model demonstrating the possi- 
bility of the extension of quantum dynamics to the case 
of nonlocal-in-time interactions. As has been shown [12], 
there are no physical reasons to restrict ourselves to the 
case of local interactions where the generalized dynamical 
equation is equivalent to the Schrodinger equation: The 
dynamics corresponding to any solution of this equation 
is not at variance with the current concepts of quantum 
physics. This means that the situation where the dynam- 
ics is generated by nonlocal-in-time interaction is possi- 
ble in principle. In the present paper we show that this 
possibility is realized in the case of low energy nucleon 
dynamics, and in leading order this dynamics is governed 
by a nonlocal-in-time interaction operator that is consid- 
ered in the exactly solvable model [12,14]. We will show 
that this feature of the GQD permits the parametriza- 
tion of the NN interactions by operators that are as well 
defined as, for example, the Coulomb potential, and are 
independent of renormalization schemes. The advantages 
of the GQD is that it allows one to describe the evolution 
of nucleon systems in a consistent way by using equations 
that do not require renormalization. In the present pa- 
per this fact will be proved in leading order of the EFT 
approach. 

In Sec. II we review the principal features of the GQD. 
In Sec. Ill we consider the exactly solvable model with 
nonlocal-in-time interaction and its applications to the 
description of the NN interaction. The dynamics of a 
quantum system within the EFT is investigated in Sec. 
IV. We will consider the nucleons as spinless particles at 
energies much smaller than their mass, their internal ex- 
citation energy, and the range of their interaction, and 
will investigate nucleon dynamics within the EFT ap- 
proach developed in Ref. [15]. We show that in leading 
order this dynamics is equivalent to the dynamics of the 
model [12,14] with a nonlocal-in-time interaction opera- 
tor. The dynamical situation that arises in an EFT after 
renormalization is investigated in Sec. V. Finally, in Sec. 
VI we present some concluding remarks. 



II. GENERALIZED QUANTUM DYNAMICS 

As has been shown in Ref. [|l2| , the Schrodinger equa- 
tion is not the most general dynamical equation consis- 
tent with the current concepts of quantum theory. Let us 
consider these concepts. As is well known, the canonical 
formalism is founded on the following assumptions: (i) 
The physical state of a system is represented by a vector 
(properly by a ray) of a Hilbert space, and (ii) an ob- 
servable A is represented by a Hermitian hypermaximal 
operator a. The eigenvalues a r of a give the possible 



values of A. An eigenvector \<fr ) corresponding to the 
eigenvalue a r represents a state in which A has the value 
a r . If the system is in the state 1-0), the probability P r 
of finding the value a r for A, when a measurement is 
performed, is given by 

s 

where Py r is the projection operator on the eigenmanifold 
V r corresponding to a r , and the sum E s is taken over a 
complete orthonormal set \(pr ) (s = 1,2, ...) of V r . The 
state of the system immediately after the observation is 
described by the vector Pv r \ip). 

In canonical formalism these postulates are used to- 
gether with the assumption that the time evolution of 
a state vector is governed by the Schrodinger equation. 
However, in QFT the Schrodinger equation is only of 
formal importance because of the UV divergences. Note 
in this connection that in the Feynman approach |L| to 
quantum theory this equation is not used as a fundamen- 
tal dynamical equation. As is well known, the main pos- 
tulate, on which this approach is founded, is as follows: 
(iii) The probability of an event is the absolute square of 
a complex number called the probability amplitude. The 
joint probability amplitude of a time-ordered sequence of 
events is product of the separate probability amplitudes 
of each of these events. The probability amplitude of an 
event which can happen in several different ways is a sum 
of the probability amplitudes for each of these ways. 

The statements of assumption (iii) express the well- 
known law for the quantum- mechanical probabilities. 
Within canonical formalism this law is derived as one 
of the consequences of the theory. However, in the Feyn- 
man formulation of quantum theory this law is used as 
the main postulate of the theory. The Feynman formula- 
tion also contains, as its essential idea, the concept of a 
probability amplitude associated with a completely spec- 
ified motion or path in space-time. From assumption 
(iii) it then follows that the probability amplitude of any 
event is a sum of the probabilities that a particle has 
a completely specified path in space-time. The contribu- 
tion from a single path is postulated to be an exponential 
whose (imaginary) phase is the classical action (in units 
of h) for the path in question. The above constitutes 
the contents of the second postulate of the Feynman ap- 
proach to quantum theory. This postulate is not so fun- 
damental as assumption (iii) , which directly follows from 
the analysis of the phenomenon of quantum interference 
|]l3| . In Ref. |T^] it has been shown that the first pos- 
tulate of the Feynman approach [assumption (iii)] can 
be used in combination with the main fundamental pos- 
tulates of canonical formalism [assumptions (i) and (ii)] 
without resorting to the second Feynman postulate and 
the assumption that the dynamics of a quantum system 
is governed by the Schrodinger equation. Such a use of 
the main assumptions of quantum theory leads to a more 
general dynamical equation than the Schrodinger equa- 
tion. 
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In a general case the time evolution of a quantum sys- 
tem is described by the evolution equation 

\*(t)) = U(t,t )Mto)), 
where U(t,to) is the unitary evolution operator, 

[/+(£, t )U(t, t Q ) = U(t, t )U+(t, t Q ) = 1, (1) 



with the group property 

U(t,t')U{t',t ) = U(t,t ), U(t ,t ) 



(2) 



Here we use the interaction picture. According to as- 
sumption (iii), the probability amplitude of an event 
which can happen in several different ways is a sum of 
contributions from each alternative way. In particular, 
the amplitude (^2|£^(£, £0)1^1) can be represented as a 
sum of contributions from all alternative ways of real- 
ization of the corresponding evolution process. Dividing 
these alternatives into different classes, we can then ana- 
lyze such a probability amplitude in different ways. For 
example, subprocesses with definite instants of the be- 
ginning and end of the interaction in the system can be 
considered as such alternatives. In this way the ampli- 
tude (tp 2 \U(t, t )\ipi) can be written in the form Jl2| 



(ih\U(t,ta)\il>i) = (Hi>i) 

t rt 2 

dt 2 / dti{ih\s(t3,h)\ih), 

J to 



(3) 



where (ip2\S(t2,ti)\if>i) is the probability amplitude that 
if at time t\ the system was in the state l^i), then the 
interaction in the system will begin at time t\ and will 
end at time £2, and at this time the system will be in 
the state \ip2). Note that in general £(£2>£i) may be 
only an operator-valued generalized function of t\ and 
t 2 , since only U(t,t ) — 1 + J* t * dt<i J^ 2 dtiS(t 2 , £1) must 
be an operator on the Hilbert space. Nevertheless, it is 
convenient to call S(t 2 , £1) an "operator" , using this word 
in the generalized sense. In the case of an isolated system 
the operator S(t 2 ,ti) can be represented in the form fl^| 

§(t 2 ,t 1 ) = exp{iHat 2 )f(t2-t 1 )exp(-iHoti), (4) 

with Ho being the free Hamiltonian. 

As has been shown in Ref. , for the evolution oper- 
ator [/(£, to) given by Eq. (^) to be unitary for any times 
to and t, the operator S(t 2 , t±) must satisfy the following 
equation: 

{t2-h)S(t 2 ,h) = f dt A I dt 3 



(5) 



x(t 4 ^t 3 )S(t 2 ,U)S(t 3 ,h). 



This equation allows one to obtain the operators S(t 2 , t\) 
for any t\ and £ 2 , if the operators 5 (£21^1) corresponding 
to infinitesimal duration times r = t' 2 — of interaction 



are known. It is natural to assume that most of the con- 
tribution to the evolution operator in the limit £2 ~~ * £1 
comes from the processes associated with the fundamen- 
tal interaction in the system under study. Denoting this 
contribution by Hi nt {t 2 ,ti), we can write 

S{t 2 ,h) -» ff <ni (£ 2) £i) + o(r e ), (6) 

where r = t 2 — t\. The parameter e is determined by 
demanding that Hi nt (t 2 , £1) must be so close to the solu- 
tion of Eq. (|5|) in the limit t 2 — > t\ that this equation has 
a unique solution having the behavior ^) near the point 
t 2 = t\. Thus this operator must satisfy the condition 

(£2-£i)i2; n t(£2,£i) -> / dt 4 [ d£ 3 (£ 4 -£ 3 ) 

J tl 

x H in t (£2 , U)H mt (£3, £1) + o(r e+1 ) . (7) 

Note that the value of the parameter e depends on the 
form of the operator Hi nt (t 2 ,ti). Since S{t 2 ,ti) and 
Hi n t{t 2 ,ti) are only operator- valued distributions, the 
mathematical meaning of the conditions @ and (ffi) needs 
to be clarified. We will assume that condition (pF means 
that 



(*2 



dt-, 



dii5(£ 2 ,£i)|*i 



(*2 



dt? 



to 



dhHintfah^i) +o(r e+ ^), 



for any vectors |^i) and | \I/ 2 ) of the Hilbert space. Con- 
dition (0) has to be considered in the same sense. 

Within GQD the operator Hi nt (t 2 ,ti) plays the role 
the interaction Hamiltonian plays in the ordinary formu- 
lation of quantum theory: It generates the dynamics of a 
system. Being a generalization of the interaction Hamil- 
tonian, this operator is called the generalized interaction 
operator. If H int (t 2 ,ti) is specified, Eq. (||) allows one 
to find the operator S(t 2 ,ti). Formula (|J) can then be 
used to construct the evolution operator U(t,to) and ac- 
cordingly the state vector, 

|V(£)> = M*o)) + f dt 2 f 2 tfti5(t2,ti)Mto)>, (8) 

J to J to 

at any time t. Thus Eq. (|^) can be regarded as an equa- 
tion of motion for states of a quantum system. By using 
Eqs. (||) and (||), the evolution operator can be repre- 
sented in the form 



1 

(n 2 \U (t,t )\ni) = (n 2 \ni) + — 

Z7T 



dx 



exp[-i(z - E n . 2 )t] exp[i(z - E ni )t ] 
(z - E n2 )(z - E ni ) 

x{m\T(z)\ ni ), 



(9) 



where z = x + iy, y)0, and 



<« 2 |T( 



z)\n\) = 1 I dr e~xjp(izT)(n 2 \T(T)\ni) . (10) 
/o 
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From Eq. (||), for the evolution operator in the 
Schrodinger picture, we get 

i f°° 

U s (t,0) = — dxexp(-izt)G(z), (11) 
2?r J^oo 



where 



{n 2 \G(z)\ ni 



z-E ni (z - E n2 )(z - E ni ] 



(n 2 \T(z)\ ni ) 



(12) 



Eq. (11) is the well-known expression establishing the 
connection between the evolution operator and the Green 
operator G(z), and can be regarded as the definition of 
the operator G(z). 

The equation of motion (|J) is equivalent to the follow- 
ing equation for the T matrix |l^| : 

d{n 2 \T{z)\n 1 ) _ \ - (n 2 \T (z)\n) (n\T {z)\ ni ) 

dz -2s (z-E n y ' [L6) 

n y J 

with the boundary condition 
(n 2 \T(z)\ ni ) 

where 



\z\ — ► oo 



{n 2 \B{z)\ ni )+o{\z\-% (14) 



B(z) = i / drexp(izT)7J^](T), 
Jo 

(3 = 1 + e, and 

#int(*2 - *i) = exp(~iH t 2 )H int (t 2 ,ti) exp(i7J ii) 

is the generalized interaction operator in the Schrodinger 
picture. The solution of Eq. ( jl3| ) satisfies the equation 

{n 2 \T(zi)\ni) - {n 2 \T(z 2 )\n 1 ) 

= z 2 -zi)> — ? —7 — — . (15) 

„ (z 2 - E n )( Zl - E„) 

This equation in turn is equivalent to the following equa- 
tion for the Green operator: 



G{ Zl ) - G{z 2 ) = (z 2 - z 1 )G{z 2 )G(z 1 ). 



(16) 



This is the Hilbert identity, which in the Hamiltonian 
formalism follows from the fact that in this case the evo- 
lution operator (11) satisfies the Schrodinger equation, 
and hence the Green operator is of the form 



G(z) = (z-H)-\ 



(17) 



with H being the total Hamiltonian. At the same time, 
as has been shown in Ref.[[l2|, Eq. (||) and hence Eqs. 
( |l3| ) and (|l^) are unique consequences of the unitarity 
condition and the representation Eq. (^), expressing the 
Feynman superposition principle [assumption (iii)]. It 
should be noted that the evolution operator constructed 
by using the Schrodinger equation can be represented 



[16] in the form (||). Being written in terms of the op- 
erators S{t 2 ,t\), Eq. (pi) does not contain operators de- 
scribing interaction in the system. It is a relation for 
S(t 2 ,ti), which are the contributions to the evolution 
operator from the processes with defined instants of the 
beginning and end of the interaction in the systems. Cor- 
respondingly, Eqs. ( |l3| ) and ( [l5| ) are relations for the T 
matrix. A remarkable feature of Eq. (||) is that it works 
as a recurrence relation, and to construct the evolution 
operator it is sufficient to know the contributions to this 
operator from the processes with infinitesimal duration 
times of interaction, i.e. from the processes of a funda- 
mental interaction in the system. This makes it possible 
to use Eq. (||) as a dynamical equation. Its form does 
not depend on the specific feature of the interaction (the 
Schrodinger equation, for example, contains the interac- 
tion Hamiltonian). Since Eq. (||) must be satisfied in all 
the cases, it can be considered as the most general dy- 
namical equation consistent with the current concepts of 
quantum theory. All the needed dynamical information 
is contained in the boundary condition for this equation, 
i.e., in the generalized interaction operator Hi nt (t 2l ti). 
As has been shown in Ref. fl2"| , the dynamics governed 
by Eq. (||) is equivalent to the Hamiltonian dynamics in 
the case where the operator H int (t 2 ,ti) is of the form 



H int {t 2 ,t 1 )=-2iS(t 2 -t 1 )H I (t 1 ), 



(18) 



with Hi(ti) being the interaction Hamiltonian in the in- 
teraction picture. In this case the state vector \i/j(t)) 
given by Eq. (0) satisfies the Schrodinger equation, 



-iHjit)^))- 



d\m) 

dt 



The delta function S(t) in Eq. ( |18[ ) emphasizes that in 
this case the fundamental interaction is instantaneous. 
Thus the Schrodinger equation results from the general- 
ized equation of motion (^J) in the case where the inter- 
action generating the dynamics of a quantum system is 
instantaneous. At the same time, Eq. (^|) permits the 
generalization to the case where the interaction gener- 
ating the dynamics of a quantum system is nonlocal in 
time jl2], In a general case, the gen eralized interac- 
tion operator has the following form |l3 : 

H mt (t 2 ,h) = -2iS(t 2 -*i)fl>(ii) + H non (t 2 ,h), 

where the first term on the right-hand side of this equa- 
tion describes the instantaneous component of the inter- 
action generating the dynamics of a quantum system, 
while the term H non (t 2 ,ti) represents its nonlocal-in- 
time part. As has been shown, there is one-to-one corre- 
spondence between nonlocality of interaction and the UV 
behavior of the matrix elements of the evolution operator 
as a function of momenta of particles: The interaction op- 
erator can be nonlocal in time only in the case where this 
behavior is "bad", i.e., in a local theory it results in UV 
divergences. In Ref. |l(| it has been shown that after 
renormalization the dynamics of the three-dimensional 
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theory of a neutral scalar field interacting through a ip 4 
coupling is governed by the generalized dynamical Eq. 
(0) with a nonlocal-in-time interaction operator. This 
lets us expect that after renormalization the dynamics of 
an EFT is also governed by this equation with a nonlocal 
interaction operator. In the next section we will consider 
this problem by using a toy model of the separable NN 
interaction. 



III. MODELS WITH NONLOCAL-IN-TIME 
INTERACTIONS AND THE SHORT-RANGE NN 
INTERACTION 

Let us consider the evolution problem for two nonrel- 
ativistic particles in the center-of-mass. We denote the 
relative momentum by p and the nucleon mass by m. 
Assume that the generalized interaction operator in the 
Schrodinger picture H^ t (r) has the form 

<p 2 |^)(r)|p 1 )=^(p 2 )^(p 1 )/(r), 

where /(t) is some function of r, and the form factor 
<ys(p) has the following asymptotic behavior for |p| — * oo : 

^(p)Hpr, (Ipi-oo). (19) 

Let, for example, <p(p) be of the form 

^(p)Hpr+g(p), 

and in the limit |p| — > oo the function <?(p) satisfies the 
estimate g(p) = o(\p\~ s ), where S > a, 6 > §. In the 

separable case, (p 2 |5(t 2 , *i)|Pi) can be represented in the 
form 

(P2|5(t 2l tl)|pi) = ¥>*(P2MPl)*(*2,tl). 

Correspondingly, (p 2 |T(z)|pi) is of the form 

(p 2 |T(z)| Pl ) =v>*(p 2 Mpi)t(z), (20) 

where, as it follows from Eq. (|l3|), the function t(z) sat- 
isfies the equation 



dt{z) 



-*»(*) 



d 3 k \ip(k)\ 2 



dz " y "' J (2tt) 3 (z-E k ) 2 

with the asymptotic condition 



where 



t{z) - fi(z) + oQz\-f>), 

\z\ — > OO 



fi[z)=i I dTexp(izr)/(r), 
Jo 



(22) 



(23) 



and Ek = j^- The solution of Eq. (|2l|) with the initial 
condition t(a) = g a , where a € (— oo, 0), is 

d 3 k \<p(k)\ 2 



t{z) = g a [l + (z- a)g a 



(2tt)3 {z-E k ){a-E k ) 



(24) 



In the case a > A, the function t(z) tends to a constant 



as \z\ 



t{z) 



X. 



(25) 



Thus in this case the function fi(z) must also tend to 
A as \z\ — >oo. From this it follows that the only possible 
form of the function /(t) is 

/(r) = -2i\S(r) + f'(r), 

where the function /'(r) has no such a singularity at the 
point r = as does the delta function. In this case the 
generalized interaction operator H-^.(t) has the form 



p a |fl-W(r)| Pl ) = -2*A<5(tV*( P2 M Pi ), 



(26) 



and hence the dynamics generated by this operator is 
equivalent to the dynamics governed by the Schrodinger 
equation with the separable potential 



(27) 



(p 2 \H I \p 1 ) = \<p*(p 2 )<p(Pi)- 



Solving Eq. ( |2l|) with the boundary condition (|25|), we 
easily get the well-known expression for the T matrix in 
the separable-potential model 



(p 2 \T(z)\ Pl ) = (p*(p2)tp(pi) 



d 3 k \ip(k)\ 2 



(2-Kf 



E k 



(28) 

Ordinary quantum mechanics does not permit the ex- 
tension of the above model to the case a < |. Indeed, in 
the case of such a large-momentum behavior of the form 
factors <p(v) , the use of the interaction Hamiltonian given 
by Eq. (|27| ) leads to the UV divergences, i.e., the inte- 
gral in ( pq ) is not convergent. We will now show that 
the generalized dynamical Eq. (]|) allows one to extend 
this model to the case — ^ < a < \. Let us determine 
the class of the functions fi(z) and correspondingly the 
value of (3 for which Eq. (|2l]) has a unique solution hav- 
ing the asymptotic behavior (p2|). In the case a < h, the 
function t(z) given by (E4J) has the following behavior for 



(21) \z\- 



t(z) -> 6!(-z) Q -5 +b 2 (-z) 2a ~ 1 

\z\ — > OO 



od^l 2 "" 1 ), (29) 



where b\ = - 
bl[M(a)+g-i 

M{a) 



-47rcos(a7r)TO Q 2 and & 2 = 61 1 a| 2 Q 
with 



d 3 k \ip(k)\ 2 



-2<i 



(2tt)S 



E k 



The parameter b\ docs not depend on g a . This means 
that all solutions of Eq. ( pl| ) have the same leading term 
in Eq. (29), and only the second term distinguishes the 
different solutions of this equation. Thus, in order to 
obtain a unique solution of Eq. (pi"]), we must specify 
the first two terms in the asymptotic behavior of t(z) 
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for \z\ — > oo. From this it follows that the functions fi(z) 
must be of the form 



2q-1 



f 1 (z) = b 1 (-z) a -* +b 2 (-z) 

and (3 = 1 — 2a. Correspondingly, the functions /(r) must 
be of the form 



f(r) = a x 



a 2 T 



(30) 



with oi = — £i>xr -1 (| — a)exp[z(— j + and a 2 = 

&2r~ 1 (l — 2a) exp(— ian), where T(z) is the gamma func- 
tion. This means that in the case a < h the generalized 
interaction operator must be of the form 



p 2 |i^(T)| Pl > = ^(p 2 Mpi) 



x a X T 



a 2 T 



(31) 



and, as it follows from Eqs. (|2(J) and (24), for the T 
matrix we have 



Pa|T(z)|pi) = N(z)<p* ( P2 Mpi), 

l^(k)| 2 



(32) 



with 

N(z) = g a [ I !- ! : - )//„ 
where 



d 3 k 



(2tt)3 (z-E k ){a-E k ) 



g a = bl(b 1 \afi- a -b2-b\M{a)) \ 

It can be easily checked that N(z) can be represented in 
the following form 



N(z) 



-b 2 + h(-z)3- a + M(z)bl 



By using Eqs. (J9J) and (p2|), we can construct the evo- 
lution operator 

i f°° 

(p 2 \U(t, t )|pi) = (p 2 |pi) + — J dx 

exp[-i(z - E P2 )t] exp[z(z - E Pl )t ] 
(z ~ E P2 ){z - E P1 ) 

x^(p 2 )ip( Pl )N(z), (33) 

where z = x + iy, and y > 0. The evolution operator 
U(t, to) defined by Eq. (B3) is a unitary operator satisfy- 
ing the composition law (ph, provided that the parameter 
b 2 is real. 

In the case a = \ the generalized interaction operator 
must be of the form [17] 



(P2|^MIPl> = -^*(P2MPl) 



/oo 
dx exp(— izr) 
-oo 



ln(-z) ln 2 (-z)y' 



(34) 



TABLE I: The parameters of the interaction operator ob- 
tained by fitting the NN data, p = lMeV -1 . 



partial wave 


a 


Cl 


13 x p 


d x p 


b 2 x p'~ ia 


3 Si(np) 


0.499 


133.5 x 10^ 


433.8 


766.2 


4.207 x 


L S {np) 


0.499 


131.8 


356.3 


3.651 x 10 b 


4.202 x 10~ 5 


L S (pp) 


0.499 


320.0 


371.7 


6.763 x 10 b 


4.204 x 10 _1> 



the model for a < ^ was used for describ- 



A 2 

where 61 = 
InRef. |J 

ing the NN interaction at low energies. The motivation 
to use such a model for parametrization of the NN forces 
is the fact that due to the quark and gluon degrees of 
freedom the NN interaction must be nonlocal in time. 
However, because of the separation of scales the system 
of hadrons should be regarded as a closed system, i.e., 
the evolution operator must be unitary and satisfy the 
composition law. From this it follows that the dynamics 
of such a system is governed by Eq. ([5]) with a nonlocal- 
in-time interaction operator. Let us now parametrize the 
NN interaction by the generalized interaction operator 
of the form Eq. (31) with the following form factor: 



<p(p) = X(p) + Cigy(p), 



with 



x(p) 



- P 



2\~- 



1 1 

— - < a < -, 

2 2 



with gy(p) being the Yamaguchi form factor, 

1 



9y (p) = 



/? 2 



- p- 



Here d, c\, and (3 are some constants. As it is the gener- 
alization of the Yamaguchi model ]l|| to the case where 
the NN interaction is nonlocal in time, our model yields 
the NN phase shifts in good agreement with experiment 
(see Figs. 1-3). The parameters of the model are quoted 
in Table I. However, the main advantage of this model 
is that it allows one to investigate the effects of the re- 
tardation in the NN interaction caused by the existence 
of the quark and gluon degrees of freedom on nucleon 
dynamics. 

As we have seen, there is the one-to-one correspon- 
dence between the form of the generalized interaction 
operator and the UV behavior of the form factor <p{p). 
In the case a > i, the operator H^t (r) must necessarily 
have the form (Eq). In this case the fundamental inter- 
action is instantaneous. In the case — h < a < | [the 
restriction a > — 5 is necessary for the integral in Eq. 

(H) to be convergent] , the only possible form of H^ t (r) 
is Eq. (|3l|), and, in the case a = ^, it must be of the form 
d34| ) , and hence the interaction generating the dynamics 
of the system is nonlocal in time. Thus the interaction 
generating the dynamics can be nonlocal in time only if 
the form factors have the " bad" large-momentum behav- 
ior that within Hamiltonian dynamics gives rise to the 
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FIG. 1: Phase shifts (solid line) in the 3 Si channel for np 
scattering, compared to the experimental data (points) (Ref. 



ultraviolet divergences: 

locality o <p(p) ~ |p|~ Q , a > ~, 

nonlocality -o ip(p) ~ |p| _c \ a < ^ - 

From this it follows that the quark-gluon retardation ef- 
fects must result in the "bad" UV behavior of the matrix 
elements of the evolution operator as a function of mo- 
menta of hadrons. Note that EFT's lead to the same con- 
clusion: Within the EFT approach the quark and gluon 
degrees of freedom manifest themselves in the form of La- 
grangians consistent with the symmetries of QCD which 
gives rise to the UV divergences. Note also that EFT's 
are local theories, despite the existence of the external 
quark and gluon degrees of freedom. However, renormal- 
ization of EFT's gives rise to the fact that these theories 
become nonlocal. 
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FIG. 2: Phase shifts (solid line) in the 1 S'o channel for np 
scattering, compared to the experimental data (points) (Ref. 
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FIG. 3: Phase shifts (solid line) in the So channel for pp 
scattering, compared to the experimental data (points) (Ref. 

{§)• 



IV. THE EFT APPROACH AND THE 
NONLOCALITY IN TIME OF EFFECTIVE 
INTERACTIONS 

Let us now investigate the character of the dynamics 
generated by the two-nucleon EFT. For the sake of sim- 
plicity we will restrict ourselves to the model in which the 
nucleons are considered as spinless. Within the EFT ap- 
proach the nucleons are considered as particles (described 
by a field 'F) with three-momenta Q much smaller than 
their mass m, the mass difference A to their first excited 
state, and the pion mass. If momenta are also small com- 
pared to the range of interaction, the NN interaction can 
be approximated by a sequence of contact interactions, 
with an increasing number of derivatives. Restricting to 
parity and time-reversal theories the effective Lagrangian 
can be written as [15] 



V ° 2m 8m 3 



V 4 



2 8 



+ -(C* 2 - C^)* + *V 2 (* + *) + . . . , (35) 
8 

where C' 2n s are parameters that depend on the details of 
the dynamics of range ~ — . The particles 'F are non- 
relativistic, and evolve only forward in time. Particles 
number is also conserved. Canonical quantization leads 
to familiar Feynman rules, and the 'F propagator at four- 
momenta p is given by 



5(p°,p) 



Pi _L P 

2m ~ 8m 3 



l£ 



(36) 



The four- 1 ]/ contact interaction is given by —iv(p,p'), 
with 



v(p,p') = C + C 2 (p 2 + p' 2 ) + 2C 2 p ■ p' 



(37) 
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with p(p') the relative momentum of the incoming (out- 
going) particles. The problem can also be easily solved 
by using the time-ordered perturbation theory, since in 
this case we deal only with the particles evolving forward 
in time. Obviously this way is more convenient for con- 
structing the off-shell T matrix. 

Let us consider the two-particle system at energy E = 

.2 i 4 

m ~~ is? + • • ■ m the center-of-mass frame. The key point 
of the EFT approach is that the problem can be solved 
by expanding in the number of derivatives at the vertices 
or particle lines. In leading order one has to keep only 
the first term Cq in Eq. (|59j), and correspondingly only 

the first two terms p° and ij~ in the propagator ( |36| ) . In 
this case the two particles evolve according to the familiar 
nonrelativistic Schrodinger propagator 



G (z) 



d 3 k |k)(k| 



(2tt)3 z - ]d 



■ IE 



(38) 



In this order the Lagrangian can be rewritten in the form 

C = *+ ( id + -!-V 2 ) * - ic* * + ** + *. (39) 
y 2m J 2 

Conservation of particle number reduces the two-nucleon 
T matrix to a sum of bubble diagrams. The ultravio- 
let divergencies can all be absorbed in the renormalized 
parameter Cq R K Summing the bubbles to a geometric 



series, one gets the T matrix Q 



(p 2 |r(°>(z)| Pl ) 



L '0 



(R) 



47T 



(40) 



Note that one can obtain the same T matrixjvia the LS 
equation with the potential Vb(p2,Pi) = C, by using 
some regularization and renormalization procedures. 

Let us now consider this problem from the point of 
view of the GQD. As we have stated, the character of the 
dynamics of the system depends on the large-momentum 
behavior of the matrix elements of the evolution oper- 
ator. Obviously, in the theory under consideration this 
behavior is determined by the dependence of the vertex 
v(p,p') on momenta. In leading order we can restrict 
ourselves to the Lagrangian ( |39"| ) and correspondingly to 
the vertex v(p,p') = Cq. This means that the T matrix 
must be of the form 



(p 2 |T( )(z)| Pl )=t(z), 



(41) 



i.e., this is the separable case with the form factor ip(p) = 
1. Taking into account that in this case a = 0, from Eq. 
(32), we get 



p 2 |T(°)(z)| Pl ) =- 



6? 



b 2 

3 

ITO2 



-l 



C. 



(R) 



47T 



(42) 



with b\ — —-^q= and Cq R ^ = Thus the requirement 

that the T matrix being of the form (41) satisfies Eq. (|l3| ) 

determines it up to one arbitrary parameter Cq R ^ , and we 
get the expression (40) which in the EFT is obtained, by 
summing the bubbles diagrams, and represents the lead- 
ing order contribution to the T matrix. The generalized 
interaction operator corresponding to the solution (42) is 
of the form 



(p 2 i4>)Ipi> 



47ri / in 

exp T 



16tt 2 



(43) 



Thus the renormalization of the EFT in leading order, i.e. 
of the theory with the Lagrangian (39), gives rise to the 
dynamics governed by the generalized dynamical Eq. (|5|) 
with the nonlocal- in-time interaction operator ( (43|) . This 
operator is a particular case of the interaction operator 
(31) of our model considered in Sec. III. Note also that, 
as has been shown in [17], the T matrix obtained in Ref. 
[19], by using the dimensional regularization of a model 
with the separable potential l^(p2,Pi) = ^ t f*{V2) t f{Vi): 



where <p(p) = (d 



P 



satisfies the generalized dy- 



namical equation with the nonlocal-in-time interaction 
operator (34). 

In leading order we have proved the fact that after 
renormalization the dynamics of the EFT under consid- 
eration is governed by the generalized dynamical Eq. (|5|) 
with a nonlocal-in-time interaction operator. We hope 
that this fact can be proved in any order of the expansion 
in the number of derivatives at the vertices or particle 
lines, and in general one can expect that renormalization 
gives rise to the fact that the dynamics of a quantum 
system is governed by the generalized dynamical Eq. (|^) 
with a nonlocal-in-time interaction operator. In Ref. [16] 
this fact has been shown, by using the example of the 
three-dimensional theory of a neutral scalar field inter- 
acting through a ip 4 coupling. The above gives reason to 
suppose that such a dynamical situation takes place af- 
ter renormalization in any theory, for example, in EFT's. 
Below we will present some general arguments leading to 
this conclusion. 

Let G\(z) be the Green operator of a renormalizablc 
theory corresponding to the momentum cutoff interac- 
tion Hamiltonian Hj (t) with renormalized constants. 
For every finite cutoff A, the operator Ga(z) obviously 
satisfies the Hilbert identity (|l6) 

Ga( Zi ) - G A (z 2 ) = (z 2 - z 1 )G A (z 2 )G A (z 1 ). (44) 

At the same time, the renormalized Green operator 
G ren (z) is a limit of the consequence of the operators 
Ga(z) for A — > oo. Since Eq. (44) is satisfied for every 
A, and contains only the operator G\(z), the renormal- 
ized Green operator must also satisfy the equation 

G ren (z\) — G ren (z 2 ) — (z 2 — zi)G ren (z 2 )G ren (zi) , 

despite the fact that the renormalized Green operator 
cannot be represented in the form (fi~7j) [in the limit 



9 



A — > oo the operators Hj(t) do not converge to some 
operator acting on the Hilbert space]. Correspondingly 
the renormalized T matrix satisfies Eqs. (|l3| ) and (15), 
despite the fact that in this case the LS and Schrodinger 
equations do not follow from these equations. Here the 
advantages of the GQD manifest themselves. Within the 
GQD the dynamical Eq. (||) is derived as a consequence 
of the most general physical principles, and for Eq. ( |l5| ) 
to be satisfied, the Green operator G(z) need not be rep- 
resented in the form (17). In the GQD this operator is 
defined by Eq. where the T matrix in turn is defined 
by Eq. (|TT|) , i.e., is expressed in terms of the amplitudes 
{ip2\S{t2,t{)\'ipi) being the contributions to the evolution 
operator from the processes in which the interaction in a 
quantum system begins at time t\ and ends at time t 2 . 
Only in the case where the interaction operator is of the 
form (|Ig|), i.e., the dynamics of the system is Hamilto- 
nian, can the operator G(z) be represented in the form 
(|l7|). For every finite cutoff A the dynamics of the sys- 
tem is Hamiltonian. At the same time, in the limiting 
case A — ► oo the dynamics is governed by the generalized 
dynamical equation with a nonlocal-in-time interaction 
operator, i.e., the dynamics is not Hamiltonian. 



V. EFFECTS OF THE NONLOCALITY OF THE 
NN INTERACTION ON THE CHARACTER OF 
NUCLEON DYNAMICS 

As we have shown, after renormalization the dynam- 
ics of the theory with Lagrangian (39), i.e., the EFT in 
leading order, is governed by the generalized dynamical 
Eq. ([5]) with a nonlocal-in-time interaction operator (43) . 
This dynamics is the same as in the model considered in 
Sec. Ill with a — 0. As has been shown in Ref. H] 
this is the case when the dynamics of a quantum system 
is not Hamiltonian. In order to clarify this point let us 
consider the specific features of the evolution operator in 
the nonlocal case a < h . In the Schrodinger picture, the 
evolution operator 



(p 2 \V(t)\ Pl ) = (P2|l7.(t,0)|pi) 

of the theory with the interaction operator (31) can be 
rewritten in the form 



(p2\V(t)\ Pl ) = < P2 | Pl )expH£ P2 <) 

i r°° exp(-»zt)(p 2 |r(s)| Pl ) 

2W-OC X (z - E P2 )(z - E pi ) ' 



(45) 



where (p 2 |T(;z)|pi) is given by Eq. (p2|). Since this T 
matrix satisfies Eqs. (|13| ) and (|1_5J), the evolution op- 
erator (33) is unitary, and satisfies the composition law 
(||). Correspondingly, the operator V(t) constitutes a 
one-parameter group of unitary operators, with the group 
property 



Assume that this group has a self-adjoint infinitesimal 
generator H which in the Hamiltonian formalism is iden- 
tified with the total Hamiltonian. Then for \ip) £ T>(H) 
we have 



t 



-iH\ 



(47) 



From this and Eq. (45) it follows that 
H = H + Hi, 



with 



(p 2 |-ff/|pi> = — / dx 



(Pain*)|pi 



(z - E P2 )(z - E pi y 



(48) 



where z 



iy, and y > 0. Since (p 2 |T(z)|pi) 



is an analytic function of z, and, in the case a < i, 
tends to zero as \z\ — >oo, from Eq. (48) it follows that 
(p 2 |ii/|pi) = for any p 2 and pi, and hence H = H . 
This means that, if the infinitesimal generator of the 
group of the operator V(t) exists, then it coincides with 
the free Hamiltonian, and the evolution operator is of the 
form V(t) — exp(-iHot). Thus, since this, obviously, is 
not true, the group of the operators V(t) has no infinites- 
imal generator, and hence the dynamics is not governed 
by the Schrodinger equation. 

It should be also noted that in the case a < h, Sfa, ti) 
is not an operator on the Hilbert space. In fact, the wave 
function 



V>(p) = (pK 

P*(p)*(t2,tl) 



= <p|s(t 2 ,ti)hk) 

r] 3 k 



(49) 



V(t 1 +t 2 )=V(t 1 )V(t 2 ), 1/(0) = 1. 



(46) 



is not square integrable for any nonzero \ipi), because of 
the slow rate of decay of the form factor ip(p) as |p| — > 00. 
Correspondingly, in the case a < \ the T matrix given by 
Eq. (32) does not represent an operator on the Hilbert 
space. However, as we have stated, in general S{t2,t{) 
may be only an operator- valued generalized function such 
that the evolution operator is an operator on the Hilbert 
space. Correspondingly, the T matrix must be such that 
G(z) given by Eq. (12) is an operator on the Hilbert 
space. The T matrix and S'(t 2 ,ii) satisfy these require- 
ments not only for a > \ but also for — i < a < ^ , since 
the evolution operator given by Eq. (33) and the corre- 
sponding G(z) are operators on the Hilbert space in this 
case. At the same time, in the case a < i we go beyond 
Hamiltonian dynamics. Thus the dynamics generated by 
the EFT is not governed by the Schrodinger equation, 
and the NN interaction cannot be parametrized by an 
interaction Hamiltonian defined on the Hilbert space. On 
the other hand, this dynamical situation is not peculiar 
from the point of view of the GQD: The Schrodinger 
equation is only a particular case of the generalized dy- 
namical Eq. (|5|) where the interaction is instantaneous, 
and the above means that the NN interaction is nonlocal 
in time. It is extremely important that within GQD the 
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EFT description of low energy nucleon dynamics becomes 
as well founded as the nonrelativistic quantum mechan- 
ics describing atomic phenomena. In fact, in both cases 
the dynamics is governed by the generalized dynamical 
equation (^), and only the interaction operators are dif- 
ferent. 

Since the interaction operator (43) represents only the 
contact component of the NN interaction, in general one 
has to consider it in combination with a long-range com- 
ponent of the interaction. In this case the interaction 
operator is of the form 



(P2|H^(r)|pi> = -2<*(r)V(p a ,p 1 ) 

+ (P2|#non(l~)|pi), 



(50) 



where (p2\H non (T)\pi) represents the contact nonlocal- 
in-time component, and V(p 2 , Pi) is a potential describ- 
ing the long-range component of the NN interaction. In 
leading order the nonlocal component is given by Eq. 
(43), and, for the interaction operator, we can write 



/ i rr (s) / \ i x 47TI / 17T 

(P2TO(r)|p 1 ) = ^r- F «p - 



16^ 2 



m 3 C, 



(R) 



-2i6(r)V(p 2 , Pl ). (51) 

The dynamical Eq. (|l3|) with the interaction operator 
(51) does not require regularization and renormalization 
and is as convenient for numerical calculations as the 
LS equation. However, in some cases it is convenient to 
reduce it to integral equations. For example, for some 
potentials the solution of the dynamical Eq. ([l3|) can be 
represented (see the Appendix) in the form 

(p 2 |T(z)| Pl ) = t (z) + h(z; Pl ) + h(z;p 2 ) 

+*2(z;P2,Pi), (52) 

where t 2 (z; P 2, Pi) is a solution of the equation 



*2(2; P 2, P l) = V(p 2 ,pi) 

tfq K{z;q,p 2 ) . , 
-i 2 (z;q,pi), 



(2nf 



z-Eo 



(53) 



with 



K(z;q,p 2 )=N(z) J ±L 



d 3 k U(p 2 ,k) 



(2tt) 3 z-E k 



^(q lP2 ), (54) 



and the functions to(z) and t\(z; p) are defined as 

d 3 k t 2 (z;p,q) 



h(z;p) = N{z) 



where 



t (z) = N{z) 1 + 



N(z) 



(2tt) 3 z-E k ' 



d 3 q fi(z;q) 
(2tt) 3 z - E q 



(55) 



(56) 



47T 



Equation (53) is a generalization of the LS equation 
to the case where the interaction operator contains a 
nonlocal- in-time component. 

In general the long-range component of the NN inter- 
action consists of the meson-exchange potentials and the 
Coulomb potential in the proton-proton channel. Since 
the contact term in Eq. (51) is of the leading order it 
is natural to keep in V(p 2 ,pi) the components of the 
same order. In Weinberg's power counting the one-pion- 
exchange potential is of the leading order. Hence in this 
order the TV TV interaction operator can be expressed as 



(p 2 |ff«(r)| Pl ) 



: exp 



7712 - 



16tt 2 



i 3 C, 



(R) 



-2iS(T)V^(p 2 , Pl ), (57) 



where V^-(p2,Pi) is the conventional one-pion-exchange 
potential. Substituting this potential into Eq. (13) and 
solving it numerically, one can easily obtain the T matrix 
and hence the evolution operator. Note that the conven- 
tional way of solving the above problem is the formal use 
of the potential [7,11] 



Vo(P2,Pi) = C + K(P2,Pi)- 



(58) 



We say "formal" since the use of such a potential leads to 
UV divergencies, and the Schrodinger and LS equations 
require regularization and renormalization. On the other 
hand, as we have shown, the contact interaction, which 
in Eq. ( [58] ) is formally represented by the term C, is 
parametrized by the operator (44) [the first two terms of 
the operator (^7|)] . In this case we deal with well-defined 
interaction operators and Eq. dl3| ) , which do not re quir e 
regularization and renormalization. By using Eq. (|12|), 
one can obtain the T matrix as easily as in the case of 
the pure one-pion-exchange potential. 

Thus the formalism of the GQD allows one to formu- 
late an EFT theory of nuclear forces as a self-consistent 
theory free from UV divergences. One of the advantages 
of such a formulation is that one can investigate the gen- 
eral consequences of the theory that have significant ef- 
fects on the character of nucleon dynamics. As is well 
known, the dynamics of many nucleon systems depends 
on the off-shell properties of the two-nucleon amplitudes. 
For this reason, let us consider the effects of the non- 
locality of the NN interaction on these properties. In 
the nonlocal case, the matrix elements of the evolution 
operator as functions of momenta do not go to zero at 
infinity as fast as is required by ordinary quantum me- 
chanics, and within Hamiltonian formalism this leads to 
ultraviolet divergences. For example, in this case the 
two-nucleon amplitudes (p2|T(z)|pi) do not go to zero 
fast enough to make the Faddeev equation well behaved. 
Another consequence of nonlocality in time of the NN in- 
teraction is that for fixed momenta pi and P2 the matrix 
elements (p 2 |T(z)|p 1 ) tend to zero as \z\ oo, while, in 
the local case, they tend to (p 2 |^|Pi) m this limit. To 
illustrate this, we present in Fig. 4 the off-shell behavior 
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FIG. 4: The behavior of the off-shell amplitude f(E) = 
10 9 (p|T(z)|p) (|p| = 500 MeV) in the 3 Si channel for np 
scattering. The solid curve corresponds to the model with a 
generalized interaction operator (fell), compared to the model 
with Yamaguchi potential with parameters given in Ref. jlj| l 
(dashed line). 

of (p2|T(z)|pi) in the limit \z\ — > oo. Thus, the non- 
locality in time of the NN interaction gives rise to an 
anomalous off-shell behavior of the two-nucleon ampli- 
tudes. The off-shell properties of the amplitudes for the 
ordinary interaction operator and the operator contain- 
ing the nonlocal term are qualitatively different. This is 
true even when the two interaction operators have ap- 
proximately the same phase shifts. 



VI. SUMMARY AND DISCUSSION 

In leading order of the EFT approach we have shown 
that the effective NN interaction is nonlocal in time, and 
low energy nucleon dynamics is governed by the dynami- 
cal Eq. (1q) with the nonlocal-in-time interaction operator 
(43). Thus, in the case of quantum mechanics of nuclc- 
ons at low energies, we deal with the non-Hamiltonian 
dynamics that has been investigated in Ref. |Q. This 
dynamics is described in a natural way by the GQD that, 
as has been shown, provides a consistent parametrization 
of the NN interaction. The leading order contact compo- 
nent of this interaction is parametrized by the nonlocal- 
in-time operator (43). This operator is well defined, 
and the generalized dynamical equation with this opera- 
tor does not requires regularization and renormalization. 
The parametrization of the NN interaction that takes 
into account the long-range component is represented by 
Eq. (51). 

The fact that in leading order of the EFT approach 
low energy nucleon dynamics is described by the model 
considered in Sec. Ill allows us to explain the advantages 
of the GQD in describing nucleon dynamics in terms of 
this model. As we have shown, the generalized dynamical 
equation permits any form factor ip(p) in Eq. (20) having 
UV behavior (19) with a > — |. In the case a > |, the 



only possible form of the interaction operator is Eq. (26). 
In this case we deal with the ordinary separable-potential 
model. In the case — \<ct<\, the interaction operator 
must be of the form (31), and hence the interaction gen- 
erating the dynamics of the system is necessarily nonlocal 
in time. In this case the T matrix is given by Eq. (32). 
On the other hand, one can obtain the same T matrix in 
another way, starting with the singular potential 

i/(p',p) = a^( p 'Mp), ^( P )H P r Q , «<^(59) 

that does not make sense without renormalization, since 
it leads to UV divergencies. By solving the LS equation 
regulated in some way and performing renormalization, 
one can get expression (32) for the T matrix. However, 
in this way one cannot determine a potential describing 
the fundamental interaction in the system. The singu- 
lar potential (|5^) is only of formal importance for the 
problem under consideration. All the information con- 
tained in this potential is that the T matrix is of the form 
(p2|T(z)|pi) = Lp* (p2)<p(pi)t(z) with the same form fac- 
tor. If the group of the evolution operators given by Eq. 
(45) had an infinitesimal generator, one could identify it 
with the Hamiltonian. However, as we have shown in Sec. 
V, this group has no infinitesimal generator, and hence 
there are no potentials that could govern the dynamics of 
the system after renormalization. Moreover, the T ma- 
trix (32) does not satisfy the LS equation and has such 
properties that are at variance with the Hamiltonian for- 
malism. Thus in this case we have only a calculation rule 
that allows one to compute results for physical quanti- 
ties. On the other hand, the above problems are the cost 
of trying to describe the dynamics of the system after 
renormalization in terms of the Hamiltonian formalism, 
despite the fact that this dynamics is non-Hamiltonian. 
From the more general point of view provided by the 
GQD we see that in this case the T matrix satisfies the 
generalized dynamical Eq. ( |l5| ) with the nonlocal-in-time 
interaction operator (31), and this operator describes the 
fundamental interaction generating the dynamics of the 
system. Thus after renormalization in the theory with 
the potential ( |59| ) we have the dynamics that, according 
to the GQD, must take a place in the case — | < a < 5, 
and an example of such a theory is low energy nucleon 
dynamics in leading order where the singular potential is 
V(P2, Pi) — C . Within the GQD this theory is internally 
consistent and is as well founded as theories with ordinary 
potentials such as quantum mechanics describing atomic 
phenomena. For example, the T matrix is well defined, 
and its properties satisfy the general requirements of the 
theory. This is not only important from the point of 
view of the internal consistency of the theory. Only a 
well-founded theory provides the possibility to obtain, in 
a theoretical way, new knowledge, and to prove the cor- 
rectness of calculations performed within its framework. 
The fact that in the approach based on the GQD we have 
the well-defined dynamical equation is a great advantage 
of this approach for practical calculations. This equa- 
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tion does not require regularization and renormalization 
and is as convenient for numerical calculations as the LS 
equation. 

It is important that the use of the GQD for investi- 
gating the dynamical situation in the EFT of nuclear 
forces gives rise to a natural parametrization of the NN 
interaction, and uniquely determines the form of the 
generalized interaction operator describing the interac- 
tion of nucleons. In contrast with initial Lagrangians of 
EFT's or formal singular potentials that are produced 
by these Lagrangians, knowing the generalized opera- 
tor of the NN interaction is sufficient to compute re- 
sults for physical quantities. We have obtained that in 
leading order the contact term of the NN interaction is 
parametrized by the interaction operator (43). In the 
same order (in Weinberg's power counting) the NN in- 
teraction is parametrized by the operator (|57|). In this 
way one can construct the interaction operators in any 
order of the EFT approach. This operator can then be 
used, for example, for determining the interaction oper- 
ators parametrizing interactions of nuclei. 



APPENDIX A 

Let us consider the solution of Eq. ( |T^ ) in the case 
where the interaction operator is of the form (51). From 
Eqs. (14) and (15) it follows that this solution can be 
represented in the form 



<P 2 |T(z)|pi 



hm (p 2 |T u (z)| Pl ) 



(Al) 



where the operator T u (z) is the solution of the equation 
T u {z) = B{u) + (u - z)B(u)G (u)G Q (z)T u {z). (A2) 
Here the operator B{z) is given by 

(p 2 |S(z)| Pl ) =/iCr) + V(pa,Pi), 

with 

4tt 16tt 2 



(A3) 



fiiz) 



The solution of Eq. (A2) can be represented in the form 

(p 2 \T u (z)\p 1 )=#\z)+t^(z;p 1 ) 

+tf i Wa) + # ) (*;pi,p 2 ). (A4) 



Substituting this representation in Eq. (A2) yields the 



following equations for t^\z), t ( "'(z;p), t\ u, (z;p 2 ) and 
H (^;Pi,P2): 

t^\z) = h{u) + {u-z)h{u) 



d 3 k 



+ti">(*;k) 



(2tt) 3 (z-E k )(u-E k ) 



(A5) 



ti(z;p) = (u- z) 



d 3 k 



(#>(*) +t?°(*;k)) 



(2tt) 3 {z-E k ){u-E k ) 

xV(k,p), (A6) 



t{ u) (z;p) = (u-z)f 1 (u) 
d 3 k (*i U) ( z ;P) +4 U) (^ k >P) 



(2tt) s 



(z-E k )(u-E k ) 



(A7) 



4 U) (^P2>Pl) = V(p 2 ,Pl) + (u-z) J 

p 1 u) (z;p 2 )+4" ) (z;p 2 ,k) 



, r w r , ^(k,Pi). (A8) 
(z - E k )(u - E k ) 

It is not difficult to verify that 

t^(z;p)=^\z;p). 

By solving the above set of equations, one can obtain 
the functions (z), t\ (z,p) and t 2 u \z, p 2 , pi) that in 
turn can be used for constructing the T matrix. In fact, 
from Eqs. (Al) and (A2) it follows that the T matrix 
can be represented in the form (|52|) where the functions 
h(z), ti(z,p) and i 2 (z,p 2 ,pi) are given by 



t (z)= lim t { Q u) (z), ti(z,p)= lim tf'(z,p), 

u — — oo u — > — oo 

*2(2,P2,Pi)= lim t 2 u \z,p 2 ,pi). 

U — * — OO 



Mi 



Taking into account that 

d 3 k u — z 



(2tt) 3 {z - E k )(u - E k ) Air 
equation ( |A5| ) can be rewritten in the form 



CM 

16^ 2 



4?r 



-llT 



ma - 




d A k 



Letting u — > — oo in this equation and assuming that 
V(p 2 ,Pi) goes to zero sufficiently fast when momenta 
tend to infinity, one can ea sily get Eq. (|5^) . In the same 
way, from Eqs. ( |A6| ) and (A8) one can derive Eqs. ( |5^ ) 
and <M). 
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